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Abstract 

We present a novel mathematical approach to model noise in dynamical systems. We do so by considering dynamics 
of a chain of diffusively coupled Nagumo cells affected by noise. We show that the noise in transmembrane current 
can be effectively modelled as fluctuations in electric characteristics of the membrane. The proposed approach to 
model noise in a nerve fibre is different from the standard additive stochastic current perturbation (the Langevin type 
equations). 
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1. Introduction 

A typical nerve fiber is coated in myelin (the myelin sheath consists of a single Schwann cell which is wrapped 
about 100 around the nerve fiber) with spatially periodic gaps, the nodes of Ranvier. Roughly the myelin sheath 
increases the membrane resistance by a factor of about 100 and decreases the membrane capacitance by a factor of 
about 100. Typically, the width of the node of Ranvier is about 1/im, the distance between nodes (the length of myelin 
sheath) is about 1.5mm that is close to 100af, where d is the nerve fiber diameter. Transmembrane ion flow occurs 
only at the nodes of Ranvier. Diffusive coupling corresponds to axial currents between nodes and allows propagation 
of changes in the transmembrane potential (action potential) in the spatial variable. Nerve fibers behave as intrinsic 
spatially discrete systems. The biological reason for such discrete structure: propagation of action potential along 
myelinated fibre is faster compare to that in nonmyelinated because of its saltatory propagation between nodes (speed 
in a myelinated fibre is aroundlOO m/s and in a nonmyelinated fiber is 1 H- 5 m/s). For further details of the model, 
physical parameters and equivalent electric circuit we refer to |Q]]. 

The idea that noise can play a positive role and benefit neural function is relatively new. Just 60 years ago it was 
commonly accepted that the noise is destructive to neural encoding J2,[3t]. Today it is well established that noise plays 
constructive role in the nerve system iQ-lsl- This new paradigm was initiated by research on stochastic resonance 
phenomena. It was shown that the stochastic resonance improves the transfer of information iMHoll. 

In this Letter, we accomplish two goals. First, we propose and study the deterministic scheme for modelling 
of noise in a nerve fiber. This scheme involves dynamical fluctuations of electric characteristics of the membrane 
together with their negative feedback control depending on the noise intensity. Then, to ensure ergodicity property 
of the dynamics, we combine this dynamical feedback control with a stochastic perturbation. In contrast to the 
random noise model (see Section|2]i our scheme operates with the only white noise process that indirectly affect initial 
dynamics. While we do not claim that our scheme to model noise in a nerve fibre is better than standard additive 
stochastic current perturbation (the Langevin type equations), we state our approach as different. 

We consider a lattice of diffusively coupled Nagumo cells described, in absence of noise, by the equations, 
iij = lAUi + f(uj), i € Z is a spacial index, where f(u) : M. — > M. has a bistable character, for example f(u) = 
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-ku(u - a)(u - 1), < a < 1, k > 0; Am, = U{+\ - 2m,- + is the standard 3 point discretization of the Laplacian 
(discrete Laplacian), and / > is a coefficient of the diffusive coupling. In addition we define "potential" V(u) by 
the differential equation, V'(u) = -f(u), V(Q) = 0. In these equations variable u corresponds to a transmembrane 
electric potential, k corresponds to the membrane conductance, a is the threshold potential. Besides modelling of the 
action potential propagation along a nerve fibre, this lattice system is important in another different areas of research 
EMI. 

In cases where i e / C Z and / is bounded, we consider this set with respect to boundary conditions, for example of 
the Neumann type. For what follows, it is convenient to represent Nagumo equations in the variational form. Define 
the "energy" functional, 'Vlu] = Yi\i\[\K^ u d 2 + V(m,)], where Vm,- = m,- - k,_i is the discrete gradient. Hereafter 
we accept short notations: ^- = 3,, jj^ = <5?, Jj = d,, and so on. With these definitions, rewrite the lattice of the 
diffusively coupled Nagumo equations in the gradient form, 

Hi = -dfV[u], i e Z. (1) 

It is easy to reveal that 'Vlu] is the Lyapunov functional since 'Vlu] = - zZ[n (dfV[u\) 2 < 0. Steady states of equations 
(Q} are the extrema of functional 'Vlu]. The minima and maxima of 'Vlu] correspond respectively to stable and 
unstable solutions of equation ([TJ. Suppose 'Vlu] > -oo and define for a continuous function A(u) the time averaging, 
A(u) = limj-^oo J Jq A(u(t))dt. Applying the time averaging to "VM = - Yj[h (dfV[u\) 2 {'Vlu] > -oo) we arrive at 

the equation, Yj{t\ (dfV[u\) 2 - 0. Thus the system spends almost all time at states of extrema of < V[u\. These extrema 
are solutions of the discrete lattice equation, /Am,- - d/'VM = 0, i e Z. This equation is implied to be equipped with 
some boundary conditions. For our purpose we accept the following conditions, Vm, — > as i — > ±oo. 

2. Random noise 

To model the influence of noise on deterministic system ([TJ, it is widely accepted in the literature that the noise is 
implemented in ([JJi by the additive stochastic currents, £;(0> ' e where {£,(f)), e z is the set of independent standard 
generalized Gaussian (5-correlated processes completely characterized by the first two cumulants, (£,(?)) = and 
j(t')\ = 5ij5(t - f'); < ■•■ > means averaging over all realizations of the random perturbations. The set of 
stochastic differential equations corresponding to (UJ takes the form, m,- = Zam, + /(m,) + V2D^,(f) = -dfV[u\ + 
y/2D^i(t), i e Z, where D is the noise intensity (we suppose that noise does not depend on node). It is convenient 
for what follows to represent this set as 

u t = - Adi'VM + yf2AD£i(t), i e Z, (2) 

where a reference time scale A is explicitly introduced. Rescaling time in (0, t — > A t, and taking into account 
scaling property of the white noise, we arrive at the case A — 1 . 

In system (f2| dissipative processes and random perturbations equilibrate one another. In respect of the "energy", 
'Vlu], we arrive at the stochastic differential equation (we specify this equation in the sense of Stratonovich (e.g. 



113)), 

<V[u] =-*Yu (WW) 2 + Vmi dMuMt)- (3) 

Equation (0 demonstrates in what way the noise affects the "energy" and defines rate of its (stochastic) fluctuations. 
Assume that 'Vlu] > -oo. Then after averaging over all realization of the random perturbations we arrive at the 
relation, -(2|;>(<WM) 2 ) + £>(£{/) d^VM) = 0, that does not depend on A; we assume {Vlu]) = const. The 



relation can be derived either by an elementary calculation or elegantly applying Novikov's formula [ 15]. This is an 
important relation that connects the noise intensity to configurational ensemble averages and thus can be considered 
as the definition of the noise intensity. In what follows we conjecture that the analogue formula involving the time 

averaging instead of the ensemble averaging, — (dfVlu]) 2 + £>2(i| d 2< Vlu] = 0, is valid and thus defines the noise 
intensity in the framework of deterministic dynamics. In order to deepen the conjecture and to describe dynamics of 
the deterministic fluctuations, we have to further presume the rate of dynamic fluctuations (r.d.f.) in the form, 

r.d.f . (dmu]) 2 + D d^lu], (4) 



that is instead of random perturbations (that do not present in deterministic dynamics) we need to consider dynamic 
fluctuations of an appropriate variable. Indeed, in absence of random perturbations we have to adopt an another way to 
properly perturb the system. Fluctuations in the electric characteristics of membrane are conjugate to that in electric 
current across the membrane. Thus it is reasonable to consider a certain electric characteristic of the membrane, 
supposedly RC (where R is resistance and C is capacitance), that defines a time scale and allows this characteristic to 
dynamically fluctuate. 

With random noise its intensity D is commonly considered as an independent parameter. Indeed, the Fokker- 
Planck operator corresponding to (|2} has the form, T*p = - <9, (dfV[u\p)+D £(,-) dfp. The Fokker-Planck equation 
associated with T* , d t p = T*p, allows the invariant solution, p M [u] ~ exp {-D _lr V[«]J. We prove the identity, 
T*p<x> [u] = 0, by straightforward calculation. It is known that this distribution and the corresponding probabilistic 
measure, dp ~ exp {-D~ lr V[w]J FT® duu are typically unique for dynamics ©. In a word, stochastic dynamics (f2j) is 

typically ergodic. This means that for every continuous function A, f A(u)dp = lim^oo j A(u(t))dt, almost for 
sure for all initial values m(0). The invariant measure relates to infinite time interval. Thus scaling the time variable 
does not affect the measure. The invariant (equilibrium) distribution p m demonstrates explicit dependence on the noise 
intensity, D. The only constraint on D arise when we presume a nondestructive role of the noise. Namely, in the case 
of a cubic nonlinearity of f(u), the general form of V(u) is double-well. Then the noise can induce transition from one 
well to another, it depends on D, and is expected to be a slow process. 

Now we can pose the problem: Given a probability measure dp (or an augmented measure on an extended phase 
space). It is necessary to find a dynamics such that lim^oo j A(u(t))dt — J A(u)dp for every continuous function 
A. We say this is dynamic modelling of the noise and assume dp ~ exp {-£T Ir V[i<]} n© duj as the invariant (ergodic) 
measure for this dynamics. 



3. Deterministic modelling of noise 

Now we will put together the above observations, - that the rate of feedback control of dynamic fluctuations and 
the invariant measure, depend on the noise intensity, - to derive a model of a deterministic noise of intensity D in a 
nerve fibre. The requirements are: 

• Dynamics of u depends on the external dynamic variables {e.g., A is endowed with its own equation of motion); 

• Rate of deterministic dynamic fluctuations is directly related to (01 (e.g., the rate of fluctuations is a measure of 
the influence of environment on electrical characteristics of the membrane); 

• Measure dp ~ exp {-D~ lr V[i<]} Y\(i) diti is invariant for the dynamics; 

• Dynamics is ergodic. 

In other words, we will sample the invariant measure, dp ~ exp {-D~ lr V[M]} dui, by the method proposed in 



lid Il7ll and to incorporate the noise intensity, D, into dynamics in accordance with This procedure is just 
reasonable since involves dynamical fluctuations of the membrane electrical characteristics. To correctly sample the 
invariant measure, dynamics must be ergodic. 

Consider dynamics in the extended phase space ({m,} , A, {rji}), 

lii — — AdfV[u\ + T]i, A — g(u), i)i — hi(u), i e Z; (5) 

functions g(u) and /i,(m) are to be determined. The extra dynamical variables A and 77, model the environment and thus 
they represent the noise effect on the Nagumo dynamics. 

Remark. Term 77, in the dynamical equations © is important. Indeed, assume 77, = 0. Then, at an equilibrium 
dfV[u\ = 0, the evolution comes to halt and no longer fluctuates, irrespective of the time dependence of A. For initial 
conditions with dfV[u\ + after a time variable rescaling, it is a gradient flow as defined in [18], and all phase space 
trajectories moves along paths with equilibrium points at either end. Thus dynamics is not ergodic. For a further 



discussion we refer to 1161 1121] 



ig path; 
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To determine functions g(u) and calculate, on the analogy of (0, 

4T«] = -A 2 (W[u]) 2 + J] i7,ama]. (6) 



(/') 



Respect to the second term on r.h.s. of (0 we put the following requirement to the time average, rjidfV[u\ — 0. A 
series of 77-dynamics satisfies this condition. Two principal limit cases are: fluctuations of current in different nodes 
are independent or synchronous. Correspondingly we endow variables {77,} with the following dynamical equations, 

i]i ~ dfV[u\, i e Z, and f] t ~ Y d/V[u], Vz € Z. (7) 

U) 

However, respect to the first term in r.h.s. of (0, we cannot repeat the trick and set A ~ (dfV[u\) 2 , since this 
results in no noise effect. To overcome this difficulty, we implement conjecture <(4j into /i-dynamics and (0, and 
explicitly set 

A ~ J] {dfVMf - D £ df<V[u] . (8) 
W {') 



Lemma 1. Asrame /I to be bounded variable, its dynamics is given by (0 and T]idfV\u\ — 0(e.g. one of dynamical 
equations (0,). 77?,en A d 2r V[u] = 0. 

Proof. First we multiply (0 by i and take into account equation (0. Then we apply the time averaging to the resulted 
equation. Thus we easily accomplish lemma. Indeed, 



£(d/VM) 2 -o2^M =-<V[u]+Y J mdmu]-DAY J dl f V[u]=DAY J tf ( VW- 

m in J © (i) (A 

□ 

This lemma together with the equations (0-(0 allows us to determine functions g(u) and hi(u) explicitly, 



g = ^Y,[(di<V[u]) 2 -Ddt<V[u]\, h i = -^dWMorh i = ~YdfV\u\, ieZ, (9) 



where Qx and Q n are parameters. Variables r\\ and corresponding functions /i, and not unique and dynamical equations 
can be simplified. 

To verify the requirement on the invariant measure we prove the theorem. 

Theorem 2. Assume the extended dynamics inform (0 where functions g and hi are given by (0, Qa > and Q n > 0. 
Then the augmented measure, 

dfi ~ exp{-D~ lr V[u]}exp[{-D~ 1 (^-Q A A 2 + ^ rj 2 )) J~~[ duidAdrji = p m J~ [ dUidAdrji, (10) 

(0 (0 © 

z's invariant for the extended dynamics. 

Remark. It should be noted that the 77-dynamics is not unique and so correspondingly it allows a variety of ?7-factors 
of the augmented measure, although they are still Gaussian. E.g., with the synchronous dynamical fluctuations, 
f] = — <9/V[h], we arrive at exp[\-D~ l jQ n ri 2 }dr]. However all cases can be treated analogously. 

Proof. The Liouville operator corresponding to the dynamics in the extended phase space has the form, £j*p = 
- di [(-Adj r V[u] + T]i)p] - dx[g(u)p] - d m [hi(u)p], and the Liouville equation reads d,p = £*p. Therefore, to 
prove the theorem we have to prove the identity, X*p<» = 0, that means the dynamics (0 preserves the augmented mea- 
sure ( TTOl l. A straightforward calculation of all partial derivatives that are involved in X*p«> with further simplification 
brings to the required identity, _£*Poo = 0. The theorem is proved. □ 



From the perspective of numerical simulations and further mathematical analysis, e.g. the Hamiltonian represen- 
tation of the proposed dynamics, it is important to find a first integral of motion. We accomplish this task with the 
following lemma. 

Lemma 3. Let the dynamical system Q and (O be augmented with the redundant dynamical variable £ £ = 
—A d?*V[zt]. Then I — ^[u] + \Q,\A 2 + jQij - D£ is the first integral of the augmented dynamical sys- 

tem. 



Proof. We derive / = by direct calculation. 



□ 



Remark. Since the origin of coordinates of the redundant variable £ is arbitrary, it is always possible for an arbitrary 
fixed trajectory to set 1 = 0. I is apparent control parameter in numerical simulations. Besides, / is related to p m and 



thus can be considered from a perspective of the Hamiltonian reformulation of dynamics on the level set / — [ 16]. 
However we do not consider this problem here. 

We can now ask whether the dynamics (01, (O is ergodic. There is no a definite answer to this question. Following 
II 1611 we can apply the Frobenius theorem of differential geometry [19] but this provides with a partial answer only. 
Here, in order to provide ergodicity, we adopt the method proposed in [16] and rigorously investigated in I20I . Namely, 
we add a Gaussian random noise to the /l-dynamics. In contrast to the model outlined in Section|2] where stochastic 
currents are added at each node, this approach relies on single and indirect stochastic perturbation. Experiments 1120 1 
reveal that, in context of the molecular dynamics, it results in a relatively weak perturbance effect on deterministic 
dynamics. Thus, we reformulate /l-dynamics (0 in the form, 



Ui = - AdfV[u] + m, A = g(u) -yA + fiyDQf&f), m = hfyi), i e Z, (11) 
where y > is a parameter. 

Theorem 4. Assume stochastically perturbed extended dynamics in the form (II lb where functions g and h are given 
by ([9j, £2/1 > 0, Q n > 0. Then the augmented measure ilOi is invariant for this dynamics. 

Proof. The Fokker-Planck operator corresponding to (HTl has the form, T*p = £j*p + yd a [(/I + DQr A l dj^ p] , and the 
Fokker-Planck equation reads d,p = T*p. After a series of routine calculations we arrive at T* p<x, = 0. Thus the 
stochastically perturbed dynamics ([Til preserves the augmented measure ( fTDl . □ 

Test simulations. Single cell dynamics. 

Low dimensional systems often reveal the ergodicity problem in a probability distribution dynamical sampling. 
For this reason, it is important to test the presented noise modelling method capable of generating the right statistic for 
a single Nagumo cell. We choose for this purpose f(u) = -4u(u - a)(u - 1). Simulations are performed using global 
parameters D = 0.04 and y = 1, for t = 10 6 . Figure 1 shows the probability distribution of the u variable calculated 
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Figure 1: Probability distributions of u variable (shown on background of exact analytical distribution). Densities are calculated as normalized 
sojourn distributions. Correspondingly, (a) a = 0.25 and (b) a = 0.5. 

with the dynamical equations and compared with exact analytical distribution. Their solid agreement brings a severe 
test of our approach. 



4. Conclusion 



We have presented a novel mathematical approach to model noise in dynamical systems. We do so by consid- 
ering dynamics of a chain of diffusively coupled Nagumo cells affected by noise. We have shown that the noise in 
transmembrane current can be effectively modelled as fluctuations in electric characteristics of the membrane. Test 
simulations give a solid support to the mathematical scheme. The proposed approach to model noise in a nerve fibre 
is different from the standard additive stochastic current perturbation and thus demonstrates a potential for further 
application. 
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